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light  intensity  at  the  observation  position  from  a  facet 

light  intensity  at  the  observation  point  from  subfacet  n 

bidirectional  reflection  distribution  function  (BRDF) 

unit  vector  bisecting  the  incident  and  reflection  directions 

surface  unit  normal 

variable  vector 

quadric  coefficient  matrix 

quadric  equation  coefficients 

center  position  of  the  facet 
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xo  observation  position 

^  angle  between  the  halfway  vector  H  and  the  incident  and  reflected 
angles 

6i  polar  angle  of  the  incident  light 

On  angle  between  the  halfway  vector  H  and  the  surface  normal 

Or  polar  angle  of  the  reflected  light 

S  microfacet  tilt  distribution  function 

a  specular  peak  width  parameter 

4* /I  luminous  power  on  a  facet 

(pi  azimuthal  angle  of  the  incident  light 

(pr  azimuthal  angle  of  the  reflected  light 


1.  Introduction 

Consider  an  object  that  is  illuminated  by  a  light  source,  such  as  a  laser.  If  the  object  is 
diffusely  reflecting,  the  light  reflected  toward  a  specific  observation  point  will  come  from 
all  parts  of  the  object  that  are  illuminated  by  the  light  source  and  observable  from  the  point 
of  view  at  the  observation  point.  If,  on  the  other  hand,  the  object  is  shiny,  most  of  the 
light  at  the  observation  point  will  come  from  the  areas  of  the  illuminated  object  where  the 
direction  of  the  specular  reflection  from  the  object  is  toward  the  observation  point.  The 
portions  of  the  object  where  the  specular  reflection  is  in  some  other  direction  will  contribute 
comparatively  little  light. 

The  most  general  way  to  model  the  light  reflecting  from  an  object  is  to  subdivide  its 
surface  into  a  large  number  of  small  facets  and  compute  the  reflection  from  each  facet.  If  a 
great  many  facets  are  used,  the  time  required  for  numerical  computations  will  be  large.  For 
a  given  observation  point,  the  faceting  does  not  need  to  be  highly  refined  except  near  the 
areas  of  the  surface  that  produce  specular  glints  toward  that  observation  point;  except  near 
the  specular  direction  the  reflected  light  is  a  weak  function  of  the  angle  and  coarse  faceting 
suffices  to  give  reasonable  accuracy.  However,  near  the  portions  of  the  object  that  produce 
specular  glints,  it  will  be  necessary  to  finely  tessellate  the  object  for  accurate  results.  Thus, 
by  faceting  most  of  the  surface  coarsely  and  using  fine  faceting  only  near  the  areas  of  glint, 
a  considerable  reduction  in  the  required  computation  can  be  achieved. 

If  one  is  considering  many  observation  points  or  a  time-dependent  situation  in  which 
the  light  source,  illuminated  object,  or  observers  are  moving,  then  different  portions  of  the 
object  will  be  areas  of  glint  for  different  observers  at  different  times.  To  ensure  that  all  of 
the  specular  glints  are  accurately  resolved,  it  will  be  necessary  to  finely  tessellate  much  or 
all  of  the  illuminated  object,  even  though  for  a  given  observer  at  a  given  time  only  a  small 
portion  of  the  surface  needs  to  be  finely  resolved. 

A  time-saving  alternative  is  to  use  adaptive  faceting:  model  the  surface  with  coarse  facets 
that  can  be  subdivided  to  more  accurately  model  the  surface  when  needed.  For  a  given  coarse 
facet,  the  level  of  subdivision  that  is  used  will  be  different  for  each  observation  point  and 
time.  Only  a  small  number  of  the  coarse  facets  will  need  to  be  significantly  subdivided, 
greatly  reducing  the  computation  time  required. 

This  adaptive  faceting  technique  has  been  used  in  the  High  Energy  Laser  Collateral 
Assessment  Tool  (HELCAT)  for  modeling  the  hazards  due  to  laser  light  reflected  off  laser 
targets.  The  technique  significantly  reduces  the  time  required  for  computational  analysis  of 
a  scenario.  i 
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Fig.  1.  A  triangle  can  be  subdivided  into  four  smaller  triangles. 


2.  Methodology 


For  reasons  to  be  discussed  later  in  this  section,  a  triangular  surface  mesh  has  been  used, 
though  similar  techniques  could  also  be  implemented  using  a  rectangular  mesh.  Consider 
a  given  coarse  facet  reflecting  light  from  a  light  source  to  a  given  observation  point  at  a 
specific  time.  One  would  first  determine  the  estimate  of  the  reflected  light  using  this  coarse 
faceting.  Consider  a  triangular  facet  defined  by  the  three  vertices  Vi,  V2,  and  V3.  Let  ^>.4 
be  the  luminous  power  on  a  facet  and  at  c  be  the  center  point  defined  as  the  position  that  is 
the  average  of  the  three  vertex  positions.  The  intensity  of  the  hght  at  the  observer  position 
xo  reflected  from  the  facet  is 


^Afr(.&i,4>i,dr<  </>r)C0S(e,)C0S(gr) 
(xc-xo)^ 


Eo  = 


(1) 


where  fr  is  the  bidirectional  reflectance  distribution  function  (BRDF),  dependent  on 
the  incident  polar  (0,)  and  azimuthal  ((/>,)  angles  and  reflected  polar  (dr)  and  azimuthal 
((f>r)  angles  as  measured  relative  to  the  normal  to  the  facet  and  some  chosen  direction 
in  the  facet  plane  indicating  the  material  anisotropy,  such  as  the  machining  grooves 
in  rolled  sheet  metal.  Here  the  facet  normal  is  taken  to  be  the  normal  to  the  plane 
defined  by  the  three  vertices  of  the  facet,  even  though  the  actual  surface  may  be 
curved. 

A  more  accurate  estimate  of  the  intensity  of  the  light  is  obtained  by  subdividing  this  initial 
coarse  facet  into  four  subfacets  as  shown  in  Fig.  1 .  The  intensity  can  then  be  computed  as 


4 


(2) 


where  Eo„  is  the  value  for  £0  as  computed  for  subfacet  n.  One  can  then  compute  the 
change  in  the  estimated  value  as  A£o  =  l^o  ~  E'q\.  If  this  change  is  small,  then  no 
further  refinement  is  needed.  If  it  is  significant,  then  each  of  the  subfacets  should  be  refined 
as  well.  The  process  can  be  continued  recursively  until  the  residual  error  is  acceptably  small 
(Fig.  2). 

Unless  the  surface  of  the  illuminated  object  is  flat,  the  subdivision  points  V12,  V23,  and 
V31  will  not  lie  in  the  plane  of  the  original  facet  and  the  subfacets  will  all  have  different 
normals.  For  the  subdivision  process  to  give  accurate  results,  an  interpolation  scheme  must 
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Fig.  2.  The  faceting  can  be  hierarchically  refined  to  as  many  levels  as  are  needed  for 
accurate  computations. 


be  available  to  determine  the  values  of  the  subdivision  points.  Two  different  schemes  have 
been  investigated:  quadric  surfaces  and  point-normal  triangles. 

A  quadric  surface  is  any  surface  defined  by  a  general  quadratic  equation  in  x,  y,  z: 

qxxX^  +  qyyy^  +  qzzZ^  +  2q^yXy  -f  2q^^xz  +  2qyzyz  +  2qjcX  +  2qyy  +  2q^z  -I-  90  =  0 

(3) 

or 


where 


and 


QP  = 

-0, 

(4) 

qxx 

qxy 

qxz 

qx' 

qxy 

9yy 

qyz 

<iy 

(5) 

qxz 

qyz 

qzz 

qz 

qx 

‘ly 

qz 

qo_ 

(6) 


Quadric  surfaces  can  represent  many  common  geometrical  shapes,  including  planes, 
cylinders,  spheres,  cones,  parabolas,  and  ellipsoids.  If  a  quadric  surface  is  associated  with 
a  coarse  facet,  it  is  simple  to  compute  the  subdivision  points  to  construct  the  subfacets. 

Most  three-dimensional  model  formats  used  by  solid  modeling  programs  represent  ob¬ 
jects  as  triangles  in  which  the  normals  at  the  vertex  points  are  separately  specified  and  not 
simply  defined  by  the  normals  to  the  flat  surface  defined  by  the  three  vertex  points.  These 
are  referred  to  as  point-normal  triangles  or  PN  triangles.  In  2001  Vlachos  et  al.  presented  a 
scheme  to  construct  a  cubic  Bezier  patch  with  quadratically  varying  normals  to  interpolate 
the  surface  position  and  normal  vectors  over  the  triangle  in  a  computationally  efficient 
manner.^  In  situations  in  which  the  actual  surface  shape  is  not  specified  for  a  triangular 
patch,  this  represents  an  excellent  means  of  approximating  the  surface. 

Because  a  quadric  surface  has  nine  independent  variables  and  three  vertices  with  three 
normals  provide  nine  constraint  equations,  it  should  be  possible  to  Solve  for  a  quadric 
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Fig.  3.  The  unit  vector  H  bisects  the  source  and  observer  directions.  The  half-angle  6s  is 
the  angle  between  H  and  the  surface  normal  n.  The  angle  between  H  and  the  source  and 
observer  directions  is  given  by  p. 


-♦-Sigma 
Sigma 
-■-Sigma 
Sigma 
—A -Sigma 
♦  Sigma 


0.03  Nonadaptive 
0.03  Adaptive 
0.01  Nonadaptive 
0.01  Adaptive 
0.003  Nonadaptive 
0.003  Adaptive 


#  of  Reflection  Calculations 

Fig.  4.  Number  of  reflection  calculations  required  when  using  adaptive  faceting  (solid  lines) 
versus  nonadaptive  faceting  (dashed  lines). 


surface  defined  by  a  set  of  vertex  points  with  specified  normals.  This  has  been  tried,  and 
sometimes  it  works  well.  However,  frequently  the  solution  matrix  is  singular  and  a  unique 
solution  cannot  be  found.  PN  triangles  do  not  require  a  matrix  solution  and  work  better  for 
modeling  cases  in  which  the  surface  shape  is  not  known  a  priori. 

3.  Example  Case 

A  test  case  was  simulated  for  light  reflection  from  a  cylindrical  object  using  the  HELCAT 
code.  The  object  was  tessellated  using  quadric  surface  elements,  allowing  the  surface  to 
approach  an  ideal  cylinder  as  the  surface  is  refined.  The  case  examined  was  for  a  fixed 
configuration  with  186  observation  points.  A  simple  but  physically  reasonable  microfacet 
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BRDF  was  used'; 


where 


4  cos  6i  ■  cos  6r  ’ 


(7) 


_  exp(-  tan^  6^ /2cr^) 
cos'*  9n 


(8) 


is  the  normalized  microfacet  tilt  distribution  function  and  9^  is  the  angle  between  the 
surface  normal  and  the  vector  bisecting  the  incident  and  reflected  directions  (Fig.  3). 

The  tests  varied  the  specular  peak  width  parameter,  a,  and  the  error  tolerance.  The  results 
are  shown  in  Fig.  4.  For  a  narrow  specular  peak  (a  =  0.(X)3),  the  adaptive  faceting  resulted 
in  an  approximately  20-fold  savings  in  the  number  of  reflection  calculations  required.  For 
a  broader  specular  peak  (a  =  0.03),  adaptive  faceting  gave  only  about  a  twofold  savings  in 
the  required  reflection  calculations. 


4.  Conclusions 

t 

Adaptive  faceting  can  considerably  reduce  the  number  of  reflection  calculations  required 
for  modeling  the  illumination  at  a  point  from  the  light  reflected  from  a  surface  with  a  narrow 
specular  peak.  In  the  absence  of  a  narrow  peak,  the  benefits  are  small. 

Similar  reductions  in  computational  requirements  might  be  possible  when  the  effective 
illumination  is  required  over  a  large  area  or  volume  by  using  a  similar  scheme  with  an 
adaptive  number  of  observation  points.  Hierarchical  refinement  would  allow  for  refining 
the  observer  spacing  only  as  required  for  specific  sets  of  reflective  facets  and  observation 
points. 
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